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Abstract. Recently we have elaborated a new cosmological model of gamma-ray 
burst (GRB) origin [1], which employs the dynamical evolution of central dense stellar 
clusters in the galactic nuclei. Those clusters inevitably contain a large fraction of 
compact stellar remnants (CSRs), such as neutron stars (NSs) and stellar mass black 
, holes (BHs), and close encounters between them result in radiative captures into short- 

£N| ' living binaries, with subsequent merging of the components, thereby producing GRBs 

(typically at large distances from the nucleus). 

In the present paper, we calculate the redshift distribution of the rate of GRBs 
produced by close encounters of NSs in distant galactic nuclei. To this end, the following 
' steps are undertaken: (i) we establish a connection between the parameters of the fast 

evolving central stellar clusters (i.e. those for which the time of dynamical evolution 
exceeds the age of the Universe) with masses of the forming central supermassive 
black holes (SMBHs) using a dynamical evolution model; (ii) we connect these masses 
with the inferred mass distributions of SMBHs in the galactic nuclei and the redshift 
distribution of quasars by assuming a certain 'Eddington luminosity phase' in their 
activity; (iii) we incorporate available observational data on the redshift distribution 
of quasars as well as a recently found correlation between the masses of galaxies and 
their central SMBHs. The resulting redshift distribution of the GRB rate, which 
accounts for both fast and slowly evolving galactic nuclei is consistent with that inferred 
from the BATSE data if the fraction of fast evolving galactic nuclei is in the range 
e~ 0.016-0.16. 



INTRODUCTION 

In our previous work [1], we have considered the generation of GRBs due to 
the radiative collisions of CSRs (for simplicity, taken to be solely NSs) in the 
central stellar clusters of evolved galactic nuclei. Depending on its initial radius 
and mass, the characteristic dynamical evolution time of a cluster, t e , can either 
exceed the age of the Universe, to, (we call such clusters as 'slowly evolving clusters' 
hereinafter) or be less than t {'fast evolving clusters 1 ). The GRB rate from slowly 



evolving clusters does not depend on redshift of the host galactic nuclei. Following 
an approach developed in [2], the mean rate of radiative collisions of NSs in the 
fiducial galactic nucleus with mass of a central cluster M, radius R, and stellar 
velocity dispersion v = (GM/2R) 1 / 2 is given by [1]: 

k(v\ 17/7 c J M \ 17/1 V R V 31/M events 
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where we use for normalization the likely values of M and R to fit the inferred rate 
of GRBs. The GRB rate from a unit comoving volume is h = (N c )n g , where the 
number density of galaxies n g ~ 1CT 2 Mpc -3 and (. . .) is an averaging throughout 
all types of slowly evolving galactic nuclei. 

A galactic nucleus with the gravitationally dominating SMBH of mass > M 
might belong to slowly evolving galactic nuclei as well. The corresponding GRB 
rate from such a nucleus is given by [1] 

9 /Nm\ 2 /v\ 17 / 7 c 

which differs from that given by Eq. (1) by a factor (M/M h ) < 1. 

One more possible contribution to GRB rate from a slowly evolving system is 
due to the coalescence of tight NS binaries in the galactic discs (which is a com- 
monly believed scheme of cosmological GRB origin). In this paper, we assume that 
there was no strong evolution with cosmological time of the number of tight NS 
binaries both in the galactic disks and in the slowly evolving nuclei so that all these 
contributions result in the redshift-independent GRB rate h ^ 1CT 8 yrs -1 Mpc -3 . 

Meanwhile the actual parameters of galactic nuclei vary in a wide range and some 
central clusters may, in fact, be fast evolving (t e < t ). The GRB rate from the 
fast evolving clusters depends on redshift. Below, we find the redshift distribution 
of the GRB rate from both slow and fast evolving clusters by (i) modelling the 
dynamical evolution of a fast evolving cluster in the galactic nucleus, which leads to 
the production of a SMBH; (ii) connecting the parameters of a fast evolving cluster 
with the mass of its central SMBH, and (iii) using the observed distributions in 
luminosity and redshift for quasars, as well as in mass for SMBHs in the galactic 
nuclei. 



FAST EVOLVING GALACTIC NUCLEI 

We introduce a dimensionless time variable y = to/t e = (l + ,2 e ) 3 / 2 , where t e is the 
time necessary for a SMBH formation, z e is the redshift at that instant. Evidently, 
y > 1 for clusters collapsed to the present epoch t — t (fast evolving clusters) and 
y < 1 for clusters which do not collapse by t (slowly evolving clusters). As we 
show elsewhere [3] , the dynamical evolution of a stellar cluster, in the framework of 
a homologous Fokker-Plank approximation [2], leads to the following relationship 



between the initial mass of the cluster and the mass of the central SMBH formed 
there: 
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There are serious observational indications [4] , [5] , [6] that a fraction of galaxies with 
a central SMBH might be as high as e ~ 0.1 and that there is a correlation between 
the central SMBH mass, M^, and the luminosity of the host galaxy bulge, L s , which 
in a simplified form can be expressed as ~ 1O~ 2 ' 5 (L S /L )M . This correlation, 
taken together with the Schechter luminosity function of galaxies [7], gives the mass 
distribution of SMBHs, <pi(Mh)dMh with 0i(M/ l ) oc e. Then we use the observed [8] 
distribution of quasars in absolute magnitude Mb and redshift z, ^(Mb, z)dMsdz, 
for z < 3. Finally, we assume the existence of an 'Eddington luminosity phase' in 
the evolution of quasars, when they shine with the Eddington luminosity (or with 
its fraction A < 1) during a certain finite time interval, which duration only depends 
upon A and the initial SMBH mass. In our calculations, the 'final' (i.e., by the end 
of the dynamical evolution) mass of a galactic nucleus is the initial mass of its new- 
born SMBH. The latter then experiences the 'Eddington luminosity phase', during 
which its mass grows up to the SMBH mass in the actually observed quasar. As 
a result, we find the mass distribution of active SMBHs, P(Mf,y)dMfdy, in terms 
of the variables y and mass Mf at the instant of the SMBH formation (in the limit 
Mf < 3 • 1O 6 M ). Thus, these active SMBHs arise due to the dynamical evolution 
of galactic nuclei with fast evolving central clusters. 



DISTRIBUTION OF GRB RATE IN REDSHIFT 

In fast evolving galactic nuclei, the resulting rate of GRB generation per unit of 
comoving volume is given as a function of redshift z by 

h ev (z) = JdM f J p( Mf ,y)N ci (M i (Mf),y) -^-^ , (4) 

where x = (1 + z) 1 ' 2 , £ ~ 0.6, N ci is found from Eq. (1) with M = Mi and R = E4, 
and function Mj(My) is the inverse function given by Eq. (3). Numerical integration 
over the variable Mf in Eq. (4) proceeds in the range Mf ~ 10 5 — 3 ■ 1O 6 M , 
which approximately correspond to the range of Mj ~ 10 7 — 2 • 1O 8 M . 

The total rate of GRBs from all kinds of galactic nuclei, h{z), includes also the 
contribution no — 10~ 8 Mpc~ 3 yr" 1 from all slowly evolving galactic nuclei and 
disks with t e > to, i.e. n(z) = no + n ev (z). The total GRB rate as a function of 
z can be fitted by a one-parametric function: n(z) = n (l + z) 13 , where the range 
of (3, 1.5 < (3 < 2, accounts for the statistical uncertainties [9]. This function is 
shown in Fig. 1. The consistency between the calculated and inferred GRB rates 
is achieved at e/X ~ 0.16. Since quasars shine at L ~ (0.1 — 1) Lgd (e.g., [10]), it 
implies that e ~ 0.016 — 0.16, which is consistent with available data [4]- [6]. 




FIGURE 1. The GRB rate as a function of redshift. The hatched region confines the allowed 
GRB redshift distribution from the BATSE data [9]. Three solid lines show the GRB rate from 
fast evolving stellar clusters in distant galactic nuclei (a constant contribution from the slowly 
evolving systems is added) as given by our model for the three different values of ej A ratio, where 
e is the fraction of galactic nuclei in which SMBHs are formed, and A is the average quasar 
luminosity normalized to the Eddington luminosity. For comparison, the dashed line shows the 
rate of GRBs from NS coalescences in tight binaries produced in the distant star-forming galaxies 
[11]. 
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